"Electro-flux" effect in superconducting hybrid Aharonov-Bohm rings 
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We have extended the circuit theory of Andreev conductance [Phys. Rev. Lett. 73, 1420 (1994)] 
to diffusive superconducting hybrid structures that contain an Aharonov-Bohm ring. The elec- 
trostatic potential distribution in the system is predicted to be flux-dependent with a period of 
the superconducting flux quantum $o = h/2e. When at least one tunnel barrier is present, the 
conductance of the system oscillates with the same period. 



Normal metal or semiconductor structures with super- 
conducting contacts have enjoyed an increasing amount 
of attention in recent years. Particularly devices known 
as Andreey interferometers have been in the focus of 
interest na The electrical transport in Andreev inter- 
ferometers depends on the phase difference of two con- 
nected superconductors, which is a clear manifestation 
of the coherent nature of multiple Andreev reflection.S 
In two recent publications&LJ two possible mechanisms 
were .discussed to explain the experiments of Petrashov 
et al.p One of them due to electron-electron interaction 
in the normal metal region and the other related to finite 
temperatures. Both mechanisms cause the resistance of 
such systems to be phase-dependent with a period of 2tt, 
in contrast to weak localization corrections to the resis- 
tance, which are predicted to display a it periodicityE 

In the present work we address a different mechanism 
that causes an oscillatory resistance in hybrid circuits. 
This effect does not depend on the phase difference be- 
tween to superconducting terminals but is due to the 
presence of a magnetic field. If a ring in the normal metal 
part of the structure is present, the voltage distribution 
and resistance is affected by a magnetic flux through the 
ring. Recently, several experiments along these lines have 
been performed!! To study this phenomenon in more de- 
tail, we will use a recently developed, easy-to-use cir- 
cuit theory of Andreev conductance. til With this theory 
is possible to calculate the zero-temperature conductance 
of diffusive hybrid systems, provided their size is small 
enough and the voltages applied are small compared to 
the magnitude of the superconducting gap. In this paper 
we extend the circuit theory in the of Ref. |^ to account 
for the presence of Aharonov-Bohm loops. We proceed by 
discussing a novel " electro-flux" effect which is in princi- 
ple present in every network that includes an Aharonov- 
Bohm ring, but is most pronounced in a circuit consisting 
solely of diffusive resistors. Although the conductance is 
in this case independent of the applied flux, the elec- 
trostatic potential distribution changes periodically with 
period <&o = The oscillatory flux-dependence of the 
conductance is computed for a few experimentally rele- 
vant geometries which include tunnel junctions. 

We consider a diffusive normal metal structure (with 
diffusion constant T>) connected to one or more super- 
conducting terminals. The circuit theory of Ref. || holds 



for sufficiently small systems: L -C £ or, equivalently, 
sufficiently small te mpera tures and voltages: T, V <C 
A, T>/L 2 . Here £ = y/D/T is the coherence length in the 
normal metal and A is the magnitude of the supercon- 
ducting gap. Finally we assume that all superconducting 
terminals are biased at the same voltage, which allows us 
to disregard non-stationary Josephson-like effects. 

The theory of Ref. | was derived using the non- 
cquilibrkwn Green function technique, originally due to 
Keldyshlia and further developed for superconductivity 
by Larkin and Ovchinnikov.tj The basic elements of 
the theory are the advanced and retarded Green func- 
tions, which determine the energy spectrum of the quasi- 
particles, and the Keldysh Green function, which de- 
scribes the filling of the spectrum by extra quasiparti- 
cles. At zero temperature, the retarded Green function 
G = s x a x + Sydy + s z a z , where a are Pauli matrices, can 
be represented by a real spectral vector s = (s x , fmS z ). 
Due to the normalization of the Green functional the 
spectral vector is also normalized: s 2 = 1. The bound- 
ary conditions on s are: s = (0, 0, 1) at all normal termi- 
nals and s = (cos <j>, sin cj>, 0) at all superconducting ones, 
where <fr equals the macroscopic phase of the supercon- 
ducting reservoir. It is thus possible to map the spatial 
phase distribution of an entire structure on the surface 
of a hemisphere. 

There are two different resistive elements, diffusive re- 
sistors and tunnel junctions. The induced superconduc- 
tivity in the normal metal region does not change the 
diffusive resistance but it does renormalizc the tunnel re- 
sistance. The expression for the spectral current (which 
is a vector in Pauli-matrix space) through a resistive el- 
ement are given by: 

D T Si X S 2 

Rnl = - = arccos(siS 2 ), (1) 

y/1 - (S1S 2 ) 2 

for a diffusive resistor with resistance Ru and 

R T I = si x sa, (2) 

for a tunnel junction with resistance s i an d s 2 are 
the spectral vectors on either side of the resistive element. 
The circuit-theory rules in terms of the spectral vectors 
are: 
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(i) The Andreev conductance of a system is the same 
as in normal circuit theory except for the fact that the 
tunnel conductivities are renormalized by a factor S1S2. 

(ii) In a normal terminal the spectral vector is the 
north pole of the hemisphere whereas in a superconduct- 
ing one it is located on the equator, where its longitude 
4> indicates the phase of the superconductor. 

(hi) The spectral current is perpendicular to both spec- 
tral vectors on either side of the resistive element. For 
a diffusive conductor the magnitude of the current is 
/ = Gdo and for a tunnel junction it is I = Gxsina. 
Here a = arccos(siS2) is the angle between the two spec- 
tral vectors at both ends of the element. 

(iv) The vector spectral current in all nodal points of 
the network is conserved. 

With these rules it is possible to compute the resistance 
of a variety of networks. However, if one wants to include 
an Aharonov-Bohm ring threaded by a flux $ into the 
circuit, these four rules have to be augmented. To see 
how this comes about we perform the standard gauge 
transformation on the Green function to get rid of the 
explicit vector potential dependence: 

G = exp(zx<r z ) G exp(-ix&z), (3) 

where \ = 77 an d $0 = 55 is the superconducting 
flux quantum. In terms of spectral vectors this gauge 
transformation is simply a rotation around the z-axis of 
the original vector by an angle of 2\- The rotated vector 
reduces to its original if 2\ — 2tt and thus will be periodic 
in the superconducting flux quantum $0. The spectral 
current vector is rotated likewise. 

Without gauge transformation (Q) the equation for G 
would be rather complicated. liil However, using (||), the 
equation for the transformed Green function G reduces 
to that for the original G in the absence of flux. The flux 
through the ring now appears in the boundary conditions 
on G as follows: At an arbitrary point P in the ring the 
Green function Gl in P and Gr infinitesimally to the 
right of P are related by: 

G L = cxp(ixa z ) G R exp(-ix<J z ), (4) 

Hence the spectral vector §l is rotated 2\ around the z- 
axis with respect to its 'neighbor' sr. Again the spectral 
current is rotated in the same way. Hence the four rules 
remain unaltered but now apply to the transformed s and 
a fifth rule is needed to prescribe the boundary condition 
in the ring: 

(v) Going around once in an Aharonov-Bohm ring, the 
spectral vector at the end of the loop is rotated by an 
angle 2\ around the z-axis with respect to the spectral 
vector at the beginning of the loop. The same holds for 
the spectral current. 

We now have all the necessary ingredients to calcu- 
late the conductance of the structures depicted in Fig. [j]. 
Network (a) consists of two diffusive wires that connect 
a normal and a superconducting terminal to a diffusive 




FIG. 1. (a) and (b) The networks under consideration. Rt 
is a tunnel junction and all other elements are diffusive resis- 
tors, (c) The circuits mapped onto a hemisphere. 

Aharonov-Bohm ring. Since a natural place for a tunnel 
barrier is at the N-S interface, we also included it in the 
circuit. Fig. |l|(b) shows a SQUID-like device, consisting 
of a ring with a tunnel junction in each branch that is 
connected to the reservoirs by two diffusive wires. We 
consider here a geometry with a single superconducting 
terminal only because we want study the effects caused 
by the applied flux rather than those due to Andreev 
interference. In Fig. 0(c) we have mapped the circuits 
onto a hemisphere to indicate the position of the spec- 
tral vectors. Because only one superconducting termi- 
nal is present, its macroscopic phase is arbitrary and we 
choose it to be zero. As can be seen from this picture we 
have chosen the point B in the ring as the point where the 
spectral vector and current are discontinuous, indicated 
schematically by the dashed line. 

Since the spectral vectors sb and sc are related by 
Eq. (|J) , we need only compute the positions of the points 
A, C and D, which are determined by spectral current 
conservation in the nodes (rule (iv)): 

E T / Gan . Gac . Gab\ n 

I = s A x ^s N -^- + s c — + s B — J =0, (5) 

£ 1 = s c x (s A ^ + s D ^) + I BC = 0, (6) 

El = s D x( Sc ^ + ^)=0, (7) 

where Gu = arccos(siSj)/y / l — (sisj) 2 and the current 
Ibc = f2 ( s b x sa) Cba/R^i where Vt is a matrix that 
rotates the spectral current over an angle —2\ according 
to rule (v). The physical current, however, is conserved 
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in every node because a uniform rotation of the spectral 
current leaves the physical current invariant. Note that 
the spectral current leaving point A is not equal to the 
spectral current arriving in C. This is a consequence of 
the gauge transformation we have used. Knowing the 
spectral vectors in the three points we are able to com- 
pute the resistance of the structure: 



Ri 



R2R3 
R2 + R3 



i?4 



cos a D s 



(8) 



where cosans = sd^s is the renormalization factor for 
tunnel conductivities according to rule (i). 

Let us now first turn to a discussion of what we call 
the electro- flux effect. Consider the geometry of Fig. 0(a) 
without the tunnel junction. In this case the total resis- 
tance of the network is not affected by the applied flux 
since the resistance of diffusive elements is not renormal- 
ized. However, the electrostatic potential distribution in 
the structure is still flux-dependent. To see this we look 
at the zem-temperature expression for the electrostatic 
potential:t3 

<p(x, $) = — Tr G K = C(x) cos 9(x, $), (9) 
4e 

where G K is the Keldysh component of the Green func- 
tion and C( x ) is t ne quasiparticle distribution-function 
that measures the deviation from equilibrium. El At zero 
temperature, Ci x ) is a linear function of position and its 
slope is proportional to the voltage drop across a resis- 
tive element. The factor cos9(x, <&) in (0) is just the 
z-component of the spectral vector, which at zero tem- 
perature is equal to the quasiparticle density of states. 
From the fact that Eq. (g) involves the flux-dependent 
9(x, 3?) it is obvious that also the electrostatic potential 
will depend on the flux through the ring. 

Fig. shows this electro-flux effect at different points 
in the structure. Here we have considered a structure 
with a total length of 3L and the wires have lengths NA 
= AC = CS = L. Fig. clearly shows that the electro- 
flux effect is largest in the middle of the structure and 
vanishes in the end points of the structure. This new 
effect is reminiscent of the electrostatic Aharonov-Bohm 
effect, in which the phase of an electron in a ring is influ- 
enced by an applied transverse electric field. Ej However, 
in a sense the electrostatic Aharonov-Bohm effect is just 
the opposite of the electro-flux effect because in the lat- 
ter case the electrostatic potential in the ring is modified 
by changing the phase of the quasiparticles with a mag- 
netic field. Using a SET-transistor it should in principle 
be possible to measure the local electrostatic potential 
in a given point. One could then measure the change in 
potential as a function of the applied flux. For a more 
detailed description of such an experiment see Ref. @. 

In the last part of this paper we discuss the flux- 
dependent conductance of several circuits that may be 
experimentally relevant. In Fig. we have plotted the 
conductance of three different systems as a function of the 
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FIG. 2. Electrostatic potential as a function of flux for dif- 
ferent points along the structure. Calculated for Ri — R4 = R 
and R 2 = R 3 = 2R. 



applied flux for different values of the resistances in the 
circuit. The conductance has been normalized to its zero- 
flux value. Panels (a) and (b) show the results for the sys- 
tem of Fig. |^(a). In panel (a) R\ = R 2 = Rj, = Ri = R 
and in panel (b) R2 = R3 = R4 = R and R\ = Rt- The 
different curves correspond to different values of R^/R. 
Fig. 0(c) shows the case in which the diffusive resistor R\ 
and the tunnel barrier i?x in Fig. 0(a) have been inter- 
changed and the remaining panel shows the conductance 
of the SQUID-like device of Fig. |(b) with R 4 = R 5 = R 
and i?i = 

Let us first consider the circuit of Fig. 0(a). The pan- 
els (a) and (b) show that, in this case, applying a flux 
through the ring decreases the conductance. This is eas- 
ily understood with the aid of Fig. 0(c). When a flux is 
applied, all points A, C and D are 'pulled' towards the 
north pole of the hemisphere, thus increasing the angle 
cvds between the spectral vectors sd and Ss- Eq. (0) 
then shows that this increases the resistance relative to 
the zero-flux value. 

It is also clear that an increase of the resistance of both 
the tunnel junction i?x and the interface resistance R\ 
causes a bigger effect on the conductance. This because 
in case (a) point D is much closer to the north pole of 
the structure than in case (b), where it is somewhere in 
the middle between N and S. Applying a flux will have a 
much larger effect on the renormalization factor cos cvds 
in case (b) than in case (a). Whereas the maximal re- 
duction in conductance in case (a) is less than a factor 
of 2, it is almost a factor of 30 in case (b). In the limit 
of large R\ and Rt our results agree with those obtained 
in Ref. § 

Although Fig. 0(a) and (b) might give the impression 
that the conductance always decreases when a flux is 
present, this is not generally the case. It is also possible 
to increase it, e.g. in systems with a single tunnel bar- 
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FIG. 3. Normalized conductance versus applied flux. Pan- 
els (a) and (b) correspond to the structure of Fig. ^(a), panel 
(c) to the same circuit with Ri and Rt interchanged and (d) 
to the SQUID-like device of Fig. 0(b). From small to large 
amplitude the different curves correspond to: (a) Rt/R—1, 2, 
3, 5, 10, 100. (b) R T /R=1, 2, 3, 5, 10, 20, 100. (c) R T /R=1, 
3, 5, 10, 20, 50, 100. (d) R T /R=3, 5, 10, 20, 50, 100. 



rier between the normal contact and the ring. When all 
diffusive resistors are kept constant and only the tunnel 
resistance is varied, the conductance is increased dramat- 
ically. As shown in Fig. |^(c) the maximum increase is a 
factor of 40 for a tunnel barrier that has a 100 times big- 
ger resistance than the diffusive resistors in the network. 

In the SQUID-like structure of Fig. |l|(b), the results 
are qualitatively the same as those shown in Fig. ||(b). 
Similar considerations as the ones used above show that 
the conductance reduction is largest when both the re- 
sistance R\ and the tunnel resistances in the ring are 
large. There is, however, a striking difference in shape 
of the curves. Whereas in panel [|(b) the minimum be- 
comes broader on increasing the resistances, the opposite 
is occurring in panel ||(d) where a sharp peak develops. 
The characteristic shapes of the curves displayed in Fig. || 
should be observable experimentally. 

In conclusion, we have generalized the circuit theory of 



Andreev conductance of Ref. || to networks that include 
an Aharonov-Bohm ring penetrated by a magnetic flux. 
We have given the complete set of altered circuit-theory 
rules and used them to calculate the flux-dependent resis- 
tance of several experimentally relevant structures. Un- 
der the right conditions these devices are very sensitive 
to the applied flux. We have predicted an electro-flux ef- 
fect in these circuits, which entails that the electrostatic 
potential distribution in the structure can be altered by 
varying the applied magnetic flux through the ring. It 
should be possible to observe this effect experimentally. 
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